Cold collisions between atoms in optical lattices 



O 

o 
o 

(N 

o 

Q 



^: 

a 

o 

-(— > 
cd 

c^ 
O 

(N ; 

>■ 

m 

^^ . 

o. 

O' 

o: 
o. 

O 

• I— I 

>^; 

Oh- 



X 



J. Piilo,^ K.-A. Suominen/'^''^ and K. Berg-S0rensen^ 

^Helsinki Institute of Physics, PL 9, FIN-OOOI4 Helsingin yliopisto, Finland 

^Department of Applied Physics, University of Turku, FIN-20014 Turun yliopisto, Finland 

^ 0rsted Laboratory, Universitetsparken 5, DK-2100 Copenhagen 0, Denmark 

*'Nordita, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 

We have simulated binary collisions between atoms in optical lattices during Sisyphus cooling. Our 
Monte Carlo Wave Function simulations show that the collisions selectively accelerate mainly the 
hotter atoms in the thermal ensemble, and thus affect the steady state which one would normally 
expect to reach in Sisyphus cooling without collisions. 
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Neutral atoms can be cooled and trapped in light- 
induced optical lattices Q . By controlling the laser light 
one can adjust the properties of the lattices in order to 
study e.g. the quantum nature of atomic motion in a pe- 
riodic structure , including the analogues to the behav- 
ior of electrons in periodic solid state lattices |^ . Ideas 
regarding the possibility to use optical lattices in atom 
optics and quantum computation have also emerged re- 
cently lQ-0. In experiments the trapped gas density is 
typically very low, providing for the moment at best a 
filling ratio of 10% for the near red detuned lattices [Q. 
Thus it is a priori a good approximation to ignore that 
the atoms interact with each other. In magneto-optical 
traps for cold atoms the inelastic collisions limit the num- 
bers and temperatures achievable for the atomic gas as 
densities increase to about 10^^ atoms/cm'^ |^. By us- 
ing Bose-Einstein condensates or combining lattices and 
other types of optical traps it is becoming possible to ob- 
tain filling ratios close to unity and even higher [pi. We 
have considered mainly the case where the filling ratio is 
about 25%, but our results can be qualitatively interpo- 
lated for smaller ratios. Applications such as quantum 
computing require atoms to interact in order to perform 
quantum logical operations ||g,g] . 

Controlled interaction studies in optical lattices could 
be performed e.g. by superimposing two optical lattices, 
which can be moved in respect to each other [g]. This, 
however, does not answer the question of what happens 
in a basic lattice configuration when the filling ratio in- 
creases, especially when inelastic collisions interfere with 
the cooling process and localization of atoms at lattice 
sites. For low densities the atom cloud reaches a thermal 
equilibrium state, and based on the studies in magneto- 
optical traps one would expect that inelastic collisions 
increase the temperature of this equilibrium state as the 
gas density increases [10|. We have performed Monte 
Carlo Wave Function (MCWF) simulations of two atoms 
in a lattice. They show that (for the parameters of our 
study) an equilibrium is not easily obtained. Instead, the 
hotter atoms are selectively accelerated, and, especially 
in 2D lattices, are likely to leave the lattice. Thus in 



densely populated lattices Sisyphus cooling could be as- 
sisted by an evaporation process: interactions eject the 
hotter atoms whereas the remaining atoms thermalize via 
Sisyphus cooling (in contrast to the coUisional thcrmal- 
ization in evaporative cooling in magnetic traps). 

In a collision two cold atoms get close enough to form a 
long-range quasimolcculc ||^. Compared to single atoms, 
the quasimolecule interacts differently with the surround- 
ing laser light, and this interaction depends on the in- 
teratomic distance. Previously the atomic interactions 
in lattices have been modelled by assuming fixed posi- 
tions for both atoms and calculating how the atomic en- 
ergy levels are shifted by the interaction |ll|-|lj]. Such 
static models ignore the dynamical nature of the inelas- 
tic collisions. But to allow the atoms to move makes the 
problem complicated and computationally tedious. We 
present in this Rapid Communication a study of colli- 
sions in a lattice between moving atoms. Once the dy- 
namical processes are understood, they can be used as 
input for macroscopic theories. This approach leaves out 
many other aspects of the problem, such as reabsorption 
of scattered photons, which is another mechanism that 
strongly limits the densities in magneto-optical traps. 
Thus our results do not necessarily reflect the complete 
situation in optical lattices, but we believe they demon- 
strate the effect of collisions on Sisyphus cooling. A more 
complete study is simply beyond the modern computa- 
tional resources. 

The distribution of atoms in an optical lattice depends 
on the choice of laser field configuration and the atomic 
level structure. The laser field should have a spatially 
changing polarization, and the atom needs at least two 
Zeeman sublevels in the lower energy state, and a differ- 
ent angular momentum in the upper energy state. The 
interaction between the laser field and an atom gives rise 
to periodic light-induced potentials for atoms in the Zee- 
man states of their internal ground states. A single atom 
moving in such a lattice will undergo Sisyphus cooling 
because of optical pumping from one ground state to an- 
other, in a manner that favors the reduction of kinetic 
energy between the rapid optical pumping cycles p5[. 



This cooling effect takes place rapidly in lattices created 
by lasers which are tuned only a few linewidths below 
the atomic transition. After cooling, the atoms are to a 
large extent localized in these potential wells. 

We have chosen as a basis for our studies the simplest 
atomic transition for a red-detuned laser field, i.e., a sys- 
tem with a lower state angular momentum Jg — 1/2 and 
an upper state angular momentum Je = 3/2. We de- 
note the first state as the ground state \g±i/2), the index 
referring to the quantum number m for the eigenvalue 
of the z projection of the angular momentum operator, 
Jz- Similarly, we denote the second state as the excited 
state, with eigenstates |e±3/2) and |e±i/2)- The reso- 
nance frequency of the transition is ujq. In the numerical 
calculations we have used the atomic properties of Cs. 

The laser field has periodicity in one dimension, and 
consists of two linearly polarized counter-propagating 
beams, with orthogonal linear polarization and frequency 
u). For this configuration, the combined laser field is 
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where So is the amplitude and k is the wavenumber. 

When the interactions become important, the atomic 
cloud is still relatively dilute so that only two atoms at 
a time are involved, and the dipole-dipole interaction 
(DDI) dominates the process. We calculate the two-atom 
DDI potentials following the procedure described in Ap- 
pendix A of Ref. [|6| . We consider two atoms interact- 
ing with the laser field, coupled to a reservoir, namely 
the vacuum electric field. The system Hamiltonian reads 
(after rotating wave approximation) 
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where the sum over a is over the two atoms, S is the 
detuning S = lu — ujq, M = 133 a.u. is the Cs atom mass. 
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The potential V gives 



the interaction with the laser field. The strength of this 
interaction is given by the Rabi frequency il = 2d£o/h 
where d is the dipole moment of the transition. 

The system interacts with the reservoir through a 
dipole coupling between the atoms and the vacuum 
modes. As we want an expression for the DDI poten- 
tial, we concentrate our effort on the calculations leading 
to an expression like Eq. (25) of Ref. flq], which origi- 
nates from the commutator between the system density 
operator, p, and the DDI potential, Vdip. Also, we con- 
centrate on spontaneous terms, i.e., terms with vanishing 
average photon number. Let us introduce the operators 
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m=-l/2 



where CC^ are the appropriate Clebsch-Gordan coeffi- 
cients and q is the polarization label in spherical basis. 



Furthermore, we use a description in terms of a center of 
mass coordinate Z and a relative coordinate r = r2 — ri 
(with coordinate along the quantization axis z — Z2 — Z1). 
With these coordinates, the interaction potential with 
the laser field reads 
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where S+^q = Sl^^ + S^^^, and A5+,g = Sl^g - S^^g. 

In order to calculate the DDI term, we look at the 
Hamiltonian part of the damping terms in the equation 
of motion for the system density operator p. After ma- 
nipulations similar to those presented in Appendix A of 
Ref. [nw, and using arguments from Ref. [n^ to evalu- 
ate integrals of Besscl functions multiplied with principal 
value functions, we find the DDI potentials between the 
two atoms. In the following we look only at atoms on 
the axis of the laser field, i.e., a one-dimensional situa- 
tion, and in this case, the DDI potential reduces to 



T/axis_ 3 Jlcosgor 



sm qar cos qor 



{S++S- 



5, 






(qorr 
O^-o) • 



(5) 



Here, F is the atomic linewidth, qo is the resonant 
wavenumber qo — uq/c, and 
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Numerical simulation of the motion of atoms in the 
lattice field in one dimension only, using the MCWF 
method |18| , is computationally very demanding |19|] . In 
order to perform two-atom studies, which require even in 
one dimension at least two translational degrees of free- 
dom, we have fixed one atom in position, and let the 
other one move freely. This fixes the relation between 
the lattice coordinates and the relative interatomic co- 
ordinate. Thus an inelastic collision will not change the 
kinetic energy for both atoms, but we use the relative ki- 
netic energy as an estimate for the kinetic energy change 
per atom. (We express energy and momentum in recoil 
units: E^ = h^k^/2M and p^ = hk respectively). 

We have formulated the problem in the two-atom basis, 
which leads to a system of 36 internal states. In studies 
for magneto-optical traps one tends to use the molecu- 
lar frame, where the atom-atom interactions have been 
included to the molecular potential structure |2^]. How- 
ever, the quantum jump processes needed for the Monte 
Carlo method are easier to describe in the atomic basis. 
One aspect of the simulations is that we do not use the 
adiabatic elimination of the excited states |2ll , which is 
typically employed in order to simplify the equations for 



atomic motion. For simplicity we neglect Doppler cool- 
ing (as Sisyphus cooling takes us below the Doppler limit, 
we expect it to be the dominant process). In the molec- 
ular frame the system of two interacting atoms is excited 
resonantly to a molecular state with an attractive inter- 
atomic potential EQ]. This leads to the acceleration of 
the relative motion of the atoms, until the process termi- 
nates with spontaneous decay. We use these attractive 
potentials for the verbal description of the process but 
it must be emphasised that they do not directly appear 
in the two-atom basis. The kinetic energy change due 
to the attractive potentials also complicates greatly the 
numerical simulations by demanding larger momentum 
and finer spatial grid than in the single atom Sisyphus 
cooling simulations |2^ . 

We use the laser parameters 5 = — SF and fl = 
1.5F, which give a lattice modulation depth of Uq = 
bSAEr- These parameters correspond to a lattice where 
the atoms move from one lattice site to another on a 
timescale that is comparable to the timescale of a har- 
monic oscillation within one of the lattice potential wells. 
In our selected system the atomic interactions are too 
weak to really destroy the lattice, so the actual case of 
interest is the one where the atoms need to be simulta- 
neously at the same lattice site. 

In the MCWF method an approximation for the two- 
atom steady state density matrix is obtained as an en- 
semble average of different wave function histories, for 
which the spontaneous emission occurs as probabilistic 



quantum jumps |18|. These quantum jumps (both atoms 
in our case have six decay channels) occur according 
to probabilities weighted by the appropriate Clebsch- 
Gordan coefficients of the decay channels. There are 
various ways how to calculate the results by ensemble 
averaging. We take the ensemble average of single his- 
tory time averages in the steady state time domain [ p3| . 
Thus we obtain the kinetic energy per atom, and the spa- 
tial and momentum probability distributions for various 
filling ratios (po) of the lattice. 

A comparison between the number of atoms having 
gained large kinetic energy via interactions and the total 
number of interaction processes show (see Table |) that 
basically every collision produces very hot atoms in our 
chosen parameter range. This leads to evaporation in 
the optical lattice: those atoms which are able to move 
from one well to the other and which have larger kinetic 
energy than localized atoms leave the trap. A crucial 
ingredient in the interaction process increasing the ki- 
netic energy by a large amount and leading to evapo- 
ration is that a large fraction of the population has to 
enter the attractive molecular excited state during the 
interaction process. This fraction in turn depends on 
the relative velocity between the interacting atoms when 
they reach the resonance point for the attractive molec- 
ular states. The relative velocity in turn depends on the 
lattice depth. In our simulations the surroundings is still 



favorable so that the relative velocity between atoms is 
low enough to keep the excitation probability high when 
atoms approach each other and cross the molecular res- 
onance point. 

The number of attractive molecular states is five and 
the resonant excitation to these potentials takes place at 
different interatomic distances [EOJ . If the atoms do not 
get a large increase in kinetic energy at the first resonance 
they reach, there are still other resonances left. A com- 
parison with semiclassical (SC) excitation and survival 
calculations suggests that the potential which becomes 
resonant first when the atoms approach each other has 
the dominant role in the inelastic collision process. 

When calculating the steady state kinetic energy per 
atom (Table ||), we use two critical wavenumbers k^- 
Wavefunction histories which at some time point have 
gained larger total kinetic energy than given by k^ are ne- 
glected in ensemble averaging (considered lost from the 
lattice). The smallest value of kc we use [g4[ is more 
than two times larger than the semiclassical critical value 
k^'^ given in Ref. |]15[. The denser the lattice is initially, 
the larger is the number of interaction processes and the 
more effective is the evaporative cooling process. This 
can be seen in the results for kinetic energy per atom 
using kc = 40 (see Table |). The kinetic energy decreases 
when the initial density of the lattice increases. Results 
with kc = 70 include atoms that are lost from the lattice, 
and the value of the kinetic energy is slightly above the 
sparse lattice (non-interacting case) result. 

The momentum distribution in Fig. ^ shows the effect 
of the evaporative cooling process clearly. Due to the 
interactions between atoms part of the population has 
shifted to the region of large k (wings in Fig. 1^) and does 
not localize back to the lattice because the atoms are 
above the recapture range. Thus the central peak of the 
momentum distribution corresponding to atoms localized 
at lattice sites has a 13% narrower FWHM for an initially 
dense lattice compared to the non-interacting case. 

We have shown that in high-density, red-detuned (a 
few linewidths) optical lattices, atomic interactions could 
lead to the ejection of the hotter atoms from the lattice, 
or at least to a selective heating process accompanied by 
a narrowing of the central momentum distribution. This 
is because (a) atoms may move from one lattice site to 
another even in the steady state for Sisyphus cooling, and 
(b) because the molecular interaction is strong enough to 
give to each clearly interacting atom pair almost always 
enough energy to escape from the lattice. Earlier simu- 
lations for roughly the same laser (and Cs) parameters 
in magneto-optical traps indicate that the dominating 
effect is a clear broadening of the atomic momentum dis- 
tribution, i.e., radiative heating llQ]. In both situations 
high-momentum atoms are produced, but in optical lat- 
tices atoms with higher momentum are strongly favoured 
in the momentum increasing process (in both cases the 
fast atoms get involved in more close encounters than the 



slow ones, but in lattices working in our selected param- 
eter region this fact becomes enhanced). 
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TABLE I. Escaped atoms and kinetic energies. The num- 
ber of MC histories which have been neglected {Nout) in the 
ensemble averaging due to escape from the lattice and steady 
state kinetic energies per atom (< Ek >) for various filling ra- 
tios (po) of the lattice. Two different critical wavenumbers kc 
have been used. N^J^t gives an estimate for the total number 
of atom-atom interaction processes based on single atom MC 
collision rate calculations by monitoring the quantum flux at 
a mean atomic separation given by po- The total number of 
MC histories for each simulation is 128. The absolute values 
of the standard deviation for the kinetic energies are given in 
parentheses. 



[1] P. S. Jessen and I. H. Deutsch, Adv. At. Mol. Opt. Phys. 
37, 95 (1996); D. R. Meacher, Contemp. Phys. 39, 329 
(1998); S. Rolston, Phys. World 11 (10), 27 (1998); L. 
Guidoni and P. Verkerk, J. Opt. B 1, R23 (1999). 

[2] Y. Castin and J. Dalibard, Europhys. Lett. 14, 761 
(1991). 

[3] G. Birkl et ai, Phys. Rev. Lett. 75, 2823 (1995); M. 
Weidemiiller et ai, ibid 75, 4583 (1995); Q. Niu et ai, 
ibid 76, 4504 (1996); M. Ben Dahan et al, ibid 76, 4508 
(1996); S. R. Wilkinson et al, ibid 76, 4512 (1996). 

[4] C. S. Adams et al, Phys. Rep. 240, 143 (1994). 

[5] G. K. Brennen et al, Phys. Rev. Lett. 82, 1060 (1999); 
G. K. Brennen et al, Phys. Rev. A 6 1, 062309 (2000); I. 



[lo; 

[11 
[12; 

[13 

[14 
[15 

[16 

[17 

[is; 

[19 



[20 
[21 
[22; 



H. Deutsch et al, :iuant-ph/0003022 



D. Jaksch et al, Phys. Rev. Lett. 82, 1975 (1999); H.-J. 

Briegel et al., J. Mod. Opt. 47, 415 (2000). 

A. Hemmerich, Phys. Rev. A 60, 943 (1999). 

K.-A. Suominen, J. Phys. B 29, 5981 (1996); J. Weiner 

et al, Rev. Mod. Phys. 71, 1 (1999). 

D. Jaksch et al, Phys. Rev. Lett. 81, 3108 (1998); D.-L 
Choi and Q. Niu, ibid. 82, 2022 (1999); M. T. DePue et 
al, ibid S2, 2262 (1999). 

M. J. Holland et ai, Phys. Rev. Lett. 72, 2367 (1994); 
Phys. Rev. A 50, 1513 (1994). 

E. V. Goldstein et al, Phys. Rev. A 53, 2604 (1996). 

C. Boisseau and J. Vigue, Opt. Commun. 127, 251 

(1996). 

A. M. Guzman and P. Meystre, Phys. Rev. A 57, 1139 

(1998). 

C. Menotti and H. Ritsch, Phys. Rev. A 60, R2653 

(1999); Appl. Phys. B 69, 311 (1999). 

J. Dalibard and C. Cohen- Tannoudji, J. Opt. Soc. Am. 

B 6, 2023 (1989); P. J. Ungar et at, ibid. 6, 2058 (1989). 

G. Lenz and P. Meystre, Phys. Rev. A 48, 3365 (1993). 

P. R. Berman, Phys. Rev. A 55, 4466 (1997). 

J. Dahbard et at, Phys. Rev. Lett. 68, 580 (1992); K. 

M0lmer et al, J. Opt. Soc. Am. B 10, 524 (1993). 

We use 32 processors of SGI Origin 2000 machine having 

128 MIPS 7?12000 processors. The total memory taken 

by a simulation is 14 Gb. Generating a single history 

requires 6 hours of CPU time. 

P. S. Julienne and J. Vigue, Phys. Rev. A 44, 4464 

(1991). 

K. I. Petsas et al, Eur. Phys. J. D 6, 29 (1999) and 

references therein. 

Due to the large kinetic energies, the spatial step size Az 



Po{%) 




Nout 


Nout 


<Ek> 


<Ek> 






kc = 40 


fc, = 70 


kc = 'iO 


kc = 70 


25.0 


39 


38 


26 


61(6) 


110(18) 


20.0 


25 


27 


19 


69(5) 


99(12) 


16.7 


19 


19 


11 


80(6) 


103(12) 


14.3 


16 


19 


12 


80(7) 


104(12) 


no interactions 











91(8) 


91(8) 



[23] 
[24] 



has to be smaller than in single atom simulations since 
the momentum space size is Lk = 2-n / Az. We use the 
split operator Fourier method to solve the non-Hermitian 
time-dependent Schrodinger equation. The lattice den- 
sity is controlled with the position space size Lz. 
K. M0lmer and Y. Castin, Quantum Semiclass. Opt. 8, 
49 (1996). 

At the semiclassical value fc^"^ |lj| the cooling force has 
its maximum value, but Sisyphus cooling is still effec- 
tive above that point. Our values of kc gives a criterion 
for neglecting energetic histories from the ensemble (i.e., 
atoms lost from the lattice). 




FIG. 1. The steady state momentum probability distribu- 
tions for a densely populated {po = 25%) lattice and for the 
non- interacting atoms case (see text). All of the MC histories 
are included. Here 5 = —3. OF, and fi = 1.5r. 



